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I describe the early days of the photon structure functions. In particular 
I discuss the parton model result of Walsh and Zerwas (1973), leading order 
QCD calculation of Witten (1976) and next-to-leading QCD calculation of 
Bardeen and myself (1978). A very brief summary of the progress made 
from 1978 until 2005 is also given. 


1. Introduction 

I have been asked to review the early days of the photon structure func¬ 
tions and to summarize the present status of this field. I have worked 
actively on deep inelastic proton and photon structure functions pQ from 
1977 to 1981, but, after summarizing the status of perturbative QCD at 
the Photon-Lepton Symposium in Bonn in 1981, I moved to study techni- 
colour models, flavour physics, weak, CP-violating and rare decays of K-and 
B-mesons, supersymmetry, extra dimensions, little Higgs and petite unifica¬ 
tion. Returning in 2005 to photon structure functions, after being decoupled 
from this field for 24 years, was a very interesting experience. One should 
realize that in 1981 no experimental data on photon structure function 
were available, although the theory had reached already a rather advanced 
stage. I felt like a person returning from a long journey to the earth. My 
main question was whether predictions for , in which I took part mainly 
in the second half of 1978, have been confirmed by the future data. I will 
give the answer to this question at the end of this writing. 

2. Simple Parton Model Result 

The story begins in the early 1970’s, when Brodsky, Kinoshita, Terezawa 
(1971), Walsh (1971), Walsh and Zerwas (February 1973) and Kingsley 
(May 1973) analyzed the deep inelastic scattering of a highly virtual pho¬ 
ton ( Q 2 ^$> p 2 ) on a real ( p 2 ~ 0) photon target as seen in Fig. |T| If the 
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photon behaved only as a vector meson, the corresponding structure func¬ 
tion i^(x,Q 2 ) would exhibit Q 2 and Bjorken x dependences similar to the 
one of the proton structure function F% (x, Q 2 ). But as pointed out by Walsh 
and Zerwas [2], at very large Q 2 the point-like contribution to F^ (x,Q 2 ), 
represented by the box diagram in Fig.Qj should dominate over the hadronic 
component. Evaluating this box diagram they found 

F^x,Q 2 ) = F^(x)log^ + ..., (2.1) 

where 

F^x) = ^^x{l-2x + 2x 2 ) (2.2) 

with ej being quark charges. I have introduced a scale A for convenience 
and the non-logarithmic terms in EU can be found in the original paper. 



Fig. 1. The basic process and the simple Parton Model. 

The most important results in m and Tn it are: 

• F^x^Q 2 ) increases at fixed x logarithmically with Q 2 as opposed to 
the corresponding decrease (at not too small x) observed for F%(x, Q 2 ). 

• The x-dependence of F^(x,Q 2 ) at large Q 2 is fully calculable as op¬ 
posed to Fo (x, Q 2 ), where only the Q 2 dependence (scaling violations) 
can be predicted within renormalization group improved perturbation 
theory but the actual shape of F^x^Q 2 ) at a given Q 2 can only be 
found by means of non-perturbative methods or directly from the data. 

3. Master Formula for F^x^Q 2 ) in QCD 

The simple parton model result of Walsh and Zerwas has been general¬ 
ized by Witten (1976) [3] to QCD in the leading logarithmic approximation 
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(LO). The NLO corrections to Witten’s result have been calculated two 
years later by Bardeen and myself (1978) |3J but the complete NNLO con¬ 
tributions have been presented for the first time only at this conference [S], 
that is 27 years later. 

The master formula for the moments of Fd(x,Q 2 ) in QCD reads as 
follows 0: 


f 


dxx n 2 F^(x, Q 2 ) = 
4ir 


a 


(Jo a s (Q 2 


a n + bn + ^ r n\ a s(Q 2 )y + ^ ^n\ a s(Q 2 )) 


,(i) 


i =1 


%— 1 


,(3.1) 


where n=2,3,.., a is the QED structure constant, a s (Q 2 ) is the QCD struc¬ 
ture constant and (3q the famous LO /3-function coefficient for which Gross, 
Politzer and Wilczek received the Nobel-Prize last year |B]. It should be 
stressed that 

(i) 

• a n , b n and ryp can be calculated in perturbation theory and correspond 
to LO, NLO and NNLO for i = 1, respectively. 

• The coefficients hn in the last sum are uncalculable in perturbative 
QCD. They can be calculated in the vector dominance model or in 
principle by means of lattice methods. They can also be extracted 
from the data. 

• As dn > 0, for very large Q 2 the first two terms dominate, except for 
one of n = 2 for which d.f' 1 = 0 and the Q 2 independent piece from 
the last sum in eu> has to be added to & 2 - 


4. More on Witten’s Analysis 

Let us recall the basic formula for the moments of the proton structure 
function: 

f dxx n ~ 2 F%(x,Q 2 ) = c h(%,9 2 )(p\Oi\p). (4.1) 

i=NS,i/>,G ^ 

Here C l n denote the Wilson coefficients (calculable in RG improved per¬ 
turbation theory) of the local operators 0^ s , O^ and Oq (quark non-singlet, 
quark singlet, gluon) and (p\0™\p) the corresponding non-perturbative ma¬ 
trix elements between the proton states. 

As pointed out by Witten 0, in the case of F^, an additional operator, 
O™, the analog of the gluon operator Oq with the non-Abelian field strength 
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tensor G a/ g replaced by the electromagnetic tensor F a p, has to be taken into 
account. The reason is that, although the Wilson coefficients Cn are 0(a), 
the matrix elements (7|0™|7) are 0(1). Therefore, the O™ contribution 
to if 2 is of the same order in a as the contribution of quark and gluon 
operators. The latter have Wilson coefficients 0(1), but matrix elements in 
photon states O(a). 

Explicitly: 


dxx n 2 F^(x,Q 2 ) = C 'n(^-^ 2 )(7|Ori7)+Q(^2,ff 2 ,a)(7|0^|7) 

(4.2) 


i=NS,ip,G 


/T 




with < 7 |0"|7> = 1. 

Including only LO contributions Witten found 


dxx 


n 2 F^(x,Q 2 ) = C^(^j,g 2 ,a) = a 2 ^ 
h Po a s (Q 


2 , Q 2 

- d n — & ^2 ' 


(4.3) 

Thus for very large Q 2 the logarithmic behaviour in Q 2 found by Walsh 
and Zerwas in the simple parton model remains valid in QCD as expected. 
However, what Witten found is that the coefficients a n of the leading loga¬ 
rithm differ from the moments of F^(x) in (HU). The anomalous behaviour 
of F^(x,Q 2 ) in QCD can be traced back to the non-vanishing anomalous 
dimension matrix in the space of the operators 0^ s , O (J and Oq, known 
already from proton deep inelastic scattering and from the anomalous di¬ 
mension Kq, describing the mixing of O™ and Oq, under renormalization. 
Setting all these anomalous dimensions to zero and keeping only the mixing 
of Ojy S and O ^ with O”, Witten’s result reduces to the one of Walsh and 
Zerwas: the anomalous behaviour of photon structure functions originates 
in certain non-vanishing anomalous dimensions. Witten’s result has been 
rederived by Llewellyn Smith (1978) using a diagrammatic method and by 
DeWitt et al.(1979) and Brodsky et al. (1978) in the framework of the 
Altarelli-Parisi approach (DGLAP). 

In summary the most important two results of Witten’s paper are the 
following ones: 


• The coefficient of logQ 2 is calculable in QCD. That is F^XjQ 2 ) is 
calculable in QCD at large Q 2 . 


• It differs from the free parton model result by finite calculable factors. 
In particular for x approaching 1 the increase of F^(x,Q 2 ) with x, 
observed also in QCD at moderate x, turns into a decrease as shown 
in Fig. [U 
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Fig. 2. FT! in the Parton Model 0, QCD at LO (3| and NLO [1|. 


5. Our 1978 NLO Analysis 

My own work in this field has been triggered by the visit of Jonathan 
Rosner to Fermilab in the summer of 1978. Jonathan wanted to know more 
about Witten’s paper and I explained it to him in the manner presented 
above. But on 24th of September 1978 I realized that Witten’s result could 
be generalized to NLO without doing a single Feynman diagram calculation. 
One had only to translate the elements of the NLO analysis for F^x^Q 2 ) 
into those for F^x^Q 2 ). A complete NLO analysis of F%(x,Q 2 ) has been 
presented for the first time in June 1978 in a work in collaboration with 
Bardeen, Duke and Muta [7]. It included the earlier result for two-loop 
anomalous dimension matrix involving O^g, O ^ and Oq of Floratos, Ross 
and Sachrajda (1977,1978). Thus in the fall of 1978 it was indeed possible 
without too much work to generalize Witten’s result to NLO. I was truly 
delighted, when Bill Bardeen agreed to join me in our second project jj]. 

The main virtues of our work, that provided the constants b n in a in 
QCD were the following ones: 

• Having b n at hand, a meaningful determination of a^ s or A jjg [7] 
from F^(x, Q 2 ) became possible. In view of the calculability of K7(x, Q 2 ) 
at large Q 2 , the prospects for the determination of a^ IS from F^(x, Q 2 ) 
appeared at least in principle better than from F^{x, Q 2 ). We will re¬ 
turn to it at the end of this writing. 

• Improved accuracy for the prediction of the shape of F^x^Q 2 ) in 
QCD. 
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• Identification of large NLO corrections for large n or equivalently large 
x, that made the turn over in the x-dependence of F^ix. Q 2 ) stronger 
than at LO as seen in Fig El 

• First comments on the divergent behaviour of the point-like compo¬ 
nent for x —> 0 that has been analyzed by many other authors in the 
1980’s. 

In the following years our result has been confirmed by many groups, 
in particular by Duke and Owens Gluck, Grassie and Reya $J and oth¬ 
ers, but in 1991 Fontannaz and Pilon Em and indepently Gluck, Reya and 
Vogt mi spotted a small error in our translation from the mixing ( Oq , Oq) 
relevant for F^ix^Q 2 ) to the mixing (0q,0™) relevant for (x,Q 2 ). This 
was a very stupid error but fortunately without essentially any numerical 
consequence for F^x^Q 2 ). A visible impact on the gluon distribution in 
the photon has been however identified. 

6. Photon Stucture Functions 1978 — 2005 

There is certainly no space to describe here in an adequate manner the 
developments after 1978. Selected reviews can be found in E3E3E1. Let 
me then list only a few points: 

On the theoretical side: 

• The evolution equations for quark and gluon distributions in the pho¬ 
ton have been studied by Brodsky et al (1978), Gluck, Grassie and 
Reya (1983), Rossi (1983), Drees (1983) and others. See additional 
references below. 

• The issue of the singular behaviour for x —> 0 has been addressed 
in several papers by Duke and Owens, Bardeen, Gliick, Grassie and 
Reya, Rossi, Antoniadis, Grunberg and others. In particular the non- 
perturbative component has been used as a regulator. 

• As stressed in particular by Gliick, Reya and Vogt (1991) fTTj . the 
Mellin n-moment technique for the study of the Q 2 -evolution of F^x, Q 2 ) 
and of quark and gluon distributions in the photon is technically supe¬ 
rior to the Bjorken-x space technique developed earlier by Gliick and 
Reya, Rossi, Drees, Da Luz Vieira, Storrow and others. In particular 
a s -counting problems can be straightforwardly avoided. 

• A large number of parametrizations of parton distributions including 
heavy flavours have been proposed. Compilations of these parametriza¬ 
tions (more than 25 in total) can be found in [121113j . Here the group 
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of Marysia Krawczyk is among the leading groups. The most recent 
efforts in these directions can be found in na naming, where nu¬ 
merous references to earlier literature can be found. 


• The NNLO corrections have been recently completed 5j. I am told 
that they are small but we have to wait until the numerical analysis 
of these corrections has been published. 


While definitely a significant progress on the study of the implications 
of the Witten’s and our calculations has been done since 1978, the main 
progress in this field has been done by experimentalists. After the first 
measurements of F^x^Q 2 ) by the CELLO collaboration in 1983 there was 
a dramatic progress in collecting data made by PLUTO, JADE, TASSO, 
TOPAZ, AMY, DELPHI, L3, ALEPH, OPAL and TPC/2. The relevant 
references can be found in [13 m H3 UHI- As a result of these efforts 
F^x^Q 2 ) and the quark distributions q' y (x,Q 2 ) are quite well known at 
present, while the gluon distribution G' r (x,Q 2 ) is still poorly known. The 
ranges in x and Q 2 explored so far are very impressive 0.001 < x < 0.9 and 
1.9GeV 2 < Q 2 < 780 GeV 2 . 



Fig. 3. as a function of Q 2 for different x compared with experiment m 
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7. Conclusions 

Let us finally ask whether the logQ 2 and x-dependences of the point¬ 
like component predicted in 1970’s have been confirmed by the data. The 
answer is given in Fig. 0 Indeed the increase of F^(x, Q 2 ) with increasing 
x and Q 2 is clearly visible. A more detailed analysis also shows that for 
x —> 1 the turn over predicted by QCD becomes visible but more data is 
required to expose this feature clearly. 

On the other hand the present data on F^(x, Q 2 ) allow already now a 
rather precise determination of a^ s . One finds m 

af s {M z ) = 0.1183 ± 0.0050(exp) ± 0.0028(theory) (7.1) 

in a very good agreement with the world average a^ ls (M z ) = 0.1182 ± 
0.0027 The accuracy could be further improved at the NLC and a photon 
collider. 

In summary, the photon structure functions predicted by QCD agree 
well with the data, even if further studies are clearly desirable. It became a 
mature field and definitely there is more to come in the next 25 years. As 
I am now returning back to flavour physics, weak decays and CP violation 
I wish all explorers of the physics of photon structure functions good luck! 
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